Size driven phase transitions in pinned vortex systems 



(N 
O 
O 
(N 

^■ 
(D 



o 
o 

u 
Oh 



X3 

o 



> 
in 

vn 

On 
O 
(N 
O 

> 



O 

o 



X 
S3 



P. S. Cornaglia and M. F. Laguna 
Centra Atomico Bariloche, 8400 San Carlos de Bariloche, Rio Negro, Argentina 

(Dated: February 1, 2008) 

We model a tridimensional vortex system in a sample with square superficial pinning in the top 
surface and obtain the ground state structures as a function of the sample thickness. Using a simple 
Frenkel-Kontorova like model and no adjustable parameters, we reproduce the experimental vortex 
configurations seen in the bottom surface and their range of stability. We find three phases with two 
transitions between them, including a continuous one from square to distorted hexagonal structure 
and a discontinuous one from distorted hexagonal to hexagonal structure. 

PACS numbers: 74.60.Ge, 74.80.-g, 74.60.Jg 



In recent years the phase diagram of vortices in 
high temperature superconductors has been the Hib- 
ject of intense experimental and theoretical studies.cl In 
Bi2Sr2CaCu208 crystals the vortex liquid transforms, 
upon cooling, for low enough magnetic fields into a 
Bragg Glass through a first-order phase transition. The 
Bragg Glass has hexagonal quasi-long range order and 
is collectively pkmed.B The presence of uncorrelated pku 
ning potentialSjS as well as periodic pinning potentials Jj'EI 
leads to the appearance of a rich variety of phases and 
transitions between them. The static and dynamic prop- 
erties of vortex systems with periodic pinnijig arrays have 
been extensively studied experimcntaUy 3l3 and, for two- 
dimensional systems, numerically.Er@lj't£ 

Superficial pinning centers with diffexent spacial dis- 
tributions can be artificially created,Ell and the mag- 
netic decoration technique allows the study of the vortex 
solid topolpgA!-pn the sample surface with a single- vortex 
resolution. PIPI In recent experiments, the structure of 
the vortex system in the presence of a superficial pin- 
ning potpatial has been studied in Bi2Sr2CaCu2 08 single 
crystals.liJ On the sample top surface, a square array of 
Fe dots is created by means of electron-beam lithography. 
These dots act as pinning centers for the vortices and are 
strong enough to determine their positions at the top 
surface in the ground state. Performing magnetic deco- 
rations on the opposite surface of the sample, with the 
magnetic field perpendicular to it, different vortex struc- 
tures as a function of the sample thickness are observed. 
For a vortex density equal to the pinning centers density 
and thick enough samples the natural hexagonal struc- 
ture is observed. This implies that the vortices deform 
themselves along the z-axis to connect the square and 
hexagonal lattices. For thin samples a square vortex ar- 
ray is found, while for intermediate sample thickness, a 
polycrystal with grains of distorted hexagonal and square 
structure is obtained. 

Despite the theoretical work done on the behavior of 
vortices in periodic potentials, to our knowledge, the 
case of the three-dimensional vortex system with two- 
dimensional pinning has not been addressed. Note that 
a direct continuous transition between square and trian- 
gular lattices, as described for example by Landau the- 
ory, cannot occur by symmetry arguments. Therefore, 



there are no obvious expectations on the evolution of the 
vortex structure with thickness. 

In this work we present a model for a three-dimensional 
vortex system in interaction with a two-dimensional pin- 
ning array at one surface of the sample and study the 
vortex configurations as a function of the sample thick- 
ness Lz- We consider the case of square superficial pin- 
ning at the first matching field (equal number of vortices 
and pinning centers) and study the strong pinning case of 
Ref. [14] , where the vortices in the ground state are pinned 
at one end. There are two competing energies, the length 
energy and the vortex-vortex interaction energy. On one 
hand, the minimum length energy is obtained for parallel 
straight vortices of length Lz- Since the vortices are fixed 
at one end to the square lattice of the top surface, this 
energy is minimized when the vortex structure is square 
at the bottom. On the other hand, the vortex-vortex 
interaction favors a triangular structure. 

For vortices tilted at small angles the voKtespiaprtex 
interaction potential can be approximated byu'taOJ 



V [u,(z), n,{z)] = {Vs K(0)] + K [u,,{Lz)]} . 



(1) 



HereMy(z) = |ui(z)-Uj(z)|, where Uj(z) = {xi{z),yi{z)) 
is the planar coordinate of the ith vortex and 



Vs (u) = 2eo 



log|-u| 



u 

T 



(2) 



This expression has a cut-off at Uc — 2.5A, where A is the 
penetration depth; u is in units of Uc and cq = (</'o/47rA)^. 
We checked that using the computationally more costly 
modified Bessel function (exact for parallel rigid vortices) 
instead of (^ gives very similar results. 

We write the length energy of the vortex i as a func- 
tion of the relative displacement of the vortex ends 
Sui = \ui{Lz) — Ui{0)\ (see Fig. |l|). Since the 5ui/Lz 
are very small, we can make a Taylor series expansion. 
The zero-th order term, ~ L^eo In k, is the energy of a 
straight vortex, where k is the Ginzburg-Landau parame- 
ter. The next non zero contribution to the length energy 
is ieo {5ui/Lz) hin. Since the results are not critically 
dependent on higher order terms, in this work we used 
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FIG. 1: Schematic section of the sample parallel to a line of 
pinning centers and to the magnetic field B. The vortex pro- 
jections, perpendicular to the plane of the figure, are drawn 
with one of their ends pinned by the superficial pinning po- 
tential. 

the following expression: 



ei{5ui) = ef)\nK^ 5ui^ + LI, (3) 

that corresponds to the length energy of a vortex tilted 
an angle arctan (Jui/Lz). We model the effect of the 
surface pinning keeping the vortex top coordinates Ui(0) 
fixed to a square lattice of lattice parameter aa . 

The total energy of the vortex system is the sum of 
the interaction energy over all the vortex pairs, and the 
length energy over all the vortices. We shift this energy 
adding a constant term, —^{Eu +E^) — Ne^, where E/\ 
and En = Lz X^i^^j [''^ii(O)] are the interaction energies 
for a triangular and square array of straight parallel vor- 
tices respectively, and N is the total number of vortices: 

r F ^ 

(4) 

This is a Frenkel-Kontorova like model in which the first 
two terms correspond to the interaction between ad- 
sorbed atoms (strain energy) and the last term models 
the effect of,the periodic substrate (elastic energy). Simi- 
lar modelsta have beeiiJjsed to study phase transitions of 
overlayers on surfaces£3 The results of this paper can be 
applied to these systems, where the decrease in the elastic 
energy -given by the substrate potential intensity- can 
be related to an increase in the number of overlayers. l3 

We generated vortex structures changing slightly the 
primitive lattice vectors of the triangular lattice in such 
a way that the vortex lattice is commensurate with both 
pinning potential and sample size. To obtain the ground 
state configurations, we performed a simulated anneal- 
ing down to zero-temperature using a hybrid algorithm. 
At each time step, we ugdate the vortex positions us- 
ing Langevin dynamics.!!!! To improve the convergence 
we allow vortex reconnections by means of a Metropo- 
lis algorithm.E3 Every five time-steps, two neighboring 
vortices {i and j) can swap their free ends Ui{Lz) and 
Uj{Lz), with a probability that depends on the energy 
difference between configurations. 

For all the ground state configurations we checked that 
the main assumption of our model remains valid, i.e. that 



5ui/Lz <C 1. Moreover, for finite and small enough tem- 
peratures this is also the case. The highest tempera- 
ture that can be studied with this model grows with de- 
creasing Lz and reaches its maximum at = 0. The 
number of vortices N is in the range [1024, 2548] and we 
use periodic boundary conditions in x and y directions. 
Here we present results for a magnetic field of 37 Gauss 
(«□ = 0.74/iTO) and k = 100. Since the disorx^r freezes 
the vortex structure at temperatures T ^ T™ J13 we have 
used the penetration length at the melting temperature 
A(T = r„i) ~ 0.8/im rather than A(T — 0), to compare 
with the experimental results. 

Hereafter we focus on the geometric structure of the 
free ends of the vortices Ui(Lz), that can be observed 
experimeritally by means of the magnetic decoration 
techniqueJU 

A free square vortex lattice is unstable under the dis- 
placement of a vortex line along the (10) or (01) lattice 
directions. Moreover, for a given displacement 5u, the 
system decreases further its energy moving neighboring 
lines of vortices in opposite directions. In an analogous 
way, when the vortices are held fixed at one end on a 
square lattice by means of a strong superficial pinning, 
they can decrease their interaction energy tilting them- 
selves. In this case, there is also an increase in the length 
of the vortices with an associated increase in the length 
energy. If the length energy increase rate, as a func- 
tion of the tilting angle, is greater than the interaction 
energy decrease rate, the square lattice with all the vor- 
tices straight is stable. This is the case for small samples 
where the length energy dominates over the interaction 
energy. 

In Fig. ||(a) we plot the ground state energies as a 
function of the sample thickness. For Lz < Ld ~ lAfim 
the elastic energy is zero and the vortices are straight in 
a square lattice. For Lz > Ld the system can decrease 
its total energy tilting the vortices. This is associated 
to an increase of the elastic energy and to a decrease 
of the strain energy. The ground state geometry changes 
at Lci from square to distorted hexagonal. The distorted 
hexagonal structure is obtained tilting lines of vortices as 
shown in the insets of Fig. ^(b). This structure is fourfold 
degenerate (the other three structures can be obtained by 
n7r/2 rotations with n integer) and can be described with 
a single parameter Su^q — Sui (in this structure the Sui 
values are the same for all vortices). For Lz < Ld the 
total energy as a function of the vortex displacement Su 
has a single minimum at Su ^ 0, as shown in the inset 
of Fig. |^(a). For Lz > Ld the Su — state becomes 
unstable, Sueq increases monotonously from zero and the 
strain energy decreases. If we keep this geometry fixed 
and let Lz ^ oo, the vortex displacements Su^q reach 
their maximum value an /3. For vortex displacements 
along the y direction, we obtain an hexagonal structure 
compressed in the x direction by a factor on/oA and 
expanded in the y direction by the inverse factor. Here 
OA is the lattice parameter of the triangular structure 
and satisfies = 2a1j/^/3. 
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FIG. 2: (a) Total {Etot), strain (Est) and elastic (Eei) ground 
state energies per vortex. Inset: Etot vs. 5u for different 
Lz ; the minimum determines Su^q (b) Equilibrium vortex dis- 
placement S a function of the sample thickness Lz . A 
sharp increase at Ld ~ lAfim is observed and corresponds 
to the transition from square to distorted hexagonal. Insets: 
square (left) and distorted hexagonal (right) structures. Open 
symbols correspond to pinning sites (top surface) and filled 
symbols to the free ends of the vortices (bottom surface). 
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FIG. 3: Typical structures for the free ends (filled symbols) 
of the vortices in the three Lz regimes. Open symbols corre- 
spond to pinning centers positions, (a) Square, (b) Distorted 
hexagonal, (c) Hexagonal, here the full vortex projection is 
also shown. Inset: detail of the structure, (d) Metastable 
state, obtained for = 4.1/im, with grains of different struc- 
ture: I) square, II) and III) distorted hexagonal in different 
orientations and IV) hexagonal. 



The square and distorted hexagonal structures are 
shown in Fig. ||(a) and Fig. ^(b) respectively. For 
the distorted hexagonal structure we show one of the 
four possible configurations. For large enough the 
distorted hexagonal structure becomes metastable. In 
fact, the minimum energy configuration for > Lc2 cor- 
responds to a rotated hexagonal structure, slightly dis- 
torted to fit the pinning geometry (inset of Fig. ^(c)), 
that for simplicity we will call hexagonal. This structure 
can be generated with the lattice vectors Si = aiOAi and 

(22 = {ai/2)aAX + (V3a2/2)aAy, where ai = 3^/^/4 ~ 
0.99 and a2 = and it is commensurate with the 

square lattice when it is rotated 45 degrees. The large 
Lz regime was also studied using a model with more in- 
ternal degrees of freedom per vortex, in which the sample 
is divided in equal thickness layers along the z direction. 
Each vortex is composed of segments with an interac- 
tion in each layer given by (0) and length energy relative 
to adjacent layers given by (H). In the range of studied 
Lz values wc find no significant difference between both 
methods. As can be seen in Fig. ||(c), where we plot the 
segment projections, the vortices deviate along the z axis 
only slightly from a straight line joining their ends at the 
sample surfaces. 

It is interesting to analyze how the symmetry break- 
ing occurs to reach the structure of Fig. ^(b) from 
high temperatures. This problem is related to the 
problem of melting in two dimensions over a square 
substrate that iia^hfien studied both theoretically and 
experimentally.c30'L3 In these systems, the substrate 
acts as an Ising perturbation for the orientational or- 
der of the overlayer, and for smooth substrates an 
Ising-like phase transition is expected as a function of 
temperatureBj For our vortex system with sample thick- 
nesses Lz ^ Lci (a "strong substrate" situation) we ob- 
serve a related behavior at finite temperatures. We find 
a temperature regime in which several grains with dif- 
ferent structures are observed: distorted hexagonal in its 
four possible orientations and square. As the tempera- 
ture decreases, the grains merge into a single crystal of 
distorted hexagonal structure (Fig. §(b)). However, if 
the vortex system is quenched, it gets trapped in a local 
energy minimum with a mixture of structures, as shown 
in Fig. ^(d). In the experiments of Ref . [14] the random 
pinning present in the saniples produces an effective slow- 
ing down in the dynamical^ and prevents the vortex sys- 
tem from reaching the ground state. For the range of 
where we find a distorted hexagonal structure, the exper- 
imental final states are a mixture of distorted hexagonal, 
square and hexagonal structures. In this regime there are 
several degenerate structures and the energy of a grain 
boundary between them vanishes at Lz = Ld- As a con- 
sequence, for finite disorder, a polycrystalline structure is 
expected with its grains decreasing in size as we approach 
Lci from above. 

At the transition from distorted hexagonal to hexag- 
onal, there is a discontinuity in the elastic and strain 
energies and in the first derivative of the total energy, as 
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FIG. 4: Ground state energies per vortex as a function of the 
sample thickness. There are three different phases separated 
by vertical bars at Lz = Ld and Lz = Lc2 ■ From left to 
right: square, distorted hexagonal and hexagonal phases. 

can be observed for — ic2 in Fig. ^. In this case, 
there is no continuous way to go from one structure to 
another through minimum energy configurations. 

As Lz is further increased we expect, in analogy with 
ID systems, a series of transitions between states with 
decreasing strain energies or eventually a commensurate- 
incommensurate transition. However, the strain energy 
difference between hexagonal lattices of different com- 
mensurations and rotation angles is very small and we 



expect the orientation and deformations to be determined 
by the random pinning, instead of the superficial pinning 
geometry. 

We showed, for the first matching field, that the system 
can be mapped to a 2D Frenkel-Kontorova model from 
the known vortex-vortex interactions and experimental 
parameters. With this simple model, we have repro- 
duced the experimental results for the vortex structure as 
a function of the sample thickness. We predict two solid- 
solid phase transitions at Lz = Ld and Lz = Lc2- The 
first one, from square to distorted hexagonal, is continu- 
ous and can be described by a single parameter Su. The 
second one is first order and corresponds to the transfor- 
mation from distorted hexagonal to hexagonal. To check 
our results we also analyzed a model with more internal 
degrees of freedom per vortex. This last model is appro- 
priate to study, in the large Lz case, how is the square to 
triangle transition along the z direction, a problem that 
will be addressed elsewhere. 

The possibility of tunning the different parameters in 
the vortex systems allows the study of Frenkel-Kontorova 
models in the presence of disorder in a controllable way. 
This systems present a rich variety of phases and transi- 
tions between them. 
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